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The possibility of achieving highly selective excitation of low metastable states of hydrogen and 
helium atoms by using short laser pulses with reasonable parameters is demonstrated theoretically. 
Interactions of atoms with the laser field are studied by solving the close-coupling equations without 
discretization. The parameters of laser pulses are calculated using different kinds of optimization 
procedures. For the excitation durations of hundreds of femtoseconds direct optimization of the 
parameters of one and two laser pulses with Gaussian envelopes is used to introduce a number of 
simple schemes of selective excitation. To treat the case of shorter excitation durations, optimal 
control theory is used and the calculated optimal fields are approximated by sequences of pulses 
with reasonable shapes. A new way to achieve selective excitation of metastable atomic states by 
using sequences of attosecond pulses is introduced. 

PACS numbers: 32.80.Qk; 32.80.Rm 



INTRODUCTION 



Nonlinear laser spectroscopy provides new possibilities 
to create and study selectively excited states of quan- 
tum systems 0, |3] ■ Development of the two-photon ex- 
citation technique Q made it possible to obtain small 
concentrations of hydrogen in the 2s metastable state. 
This was sufficient to carry out precise optical measure- 
ments of the hyperfine structure of 2s state [4] and other 
relativistic and radiative effects p. Metastable atoms 
and atomic ions also play an important role in excitation 
and charge transfer processes, even in high-temperature 
laboratory and astrophysical plasmas 0. Small concen- 
trations of metastable atoms and atomic ions can signif- 
icantly affect the radiative spectra of plasmas because of 
large cross sections of electronic capture to the excited 
states of highly-charged ions during collisions with the 
metastable atoms 0. This conclusion is based on the- 
oretical and indirect spectroscopic observations. Unfor- 
tunately, direct measurements are still very difficult due 
to the low efficiency of modern techniques for the pro- 
duction of beams of metastable atoms Q • The above ex- 
amples demonstrate the importance of the development 
of effective methods for selective excitation of metastable 
atomic states. 

Recent progress in laser techniques as well as in theo- 
retical understanding of multiphoton processes in atoms 
makes it possible to introduce novel methods for effective 
selective excitation of atomic states by using short laser 
pulses. A number of approaches to control population 
transfer between atomic and molecular states have been 
proposed so far. They can be roughly divided into two 
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groups according to the basic principles. The methods 
from the first group exploit some known physical mech- 
anisms lililiflOlilllllilllllillillil, e.g., it 

is possible to introduce some reliable and solvable model 
equations to estimate the parameters of the controlling 
laser field analytically, while the methods from the second 
group are based on the idea of direct variation of the laser 
field to maximize the desired output 19, 20, 21, 22, 23]. 

The schemes from the first group ultimately rely on 
the properties of 7r-pulses or on adiabatic properties. The 
controlling laser field is found as an accurate or approx- 
imate solution of the coupled equations within rotating 
wave approximation. These schemes, in turn, can be 
classified by the number of laser pulses involved. 

Population inversion achieved by using a single 7r-pulse 
represents the simplest case [sl . The frequency of the 
laser should be in resonance with the transition energy 
and the control is performed by changing the pulse area. 
This method can also be used to control multiphoton 
transitions. In that case, the frequency of the laser pulse 
should be a fraction of the transition energy so that the 
energy conservation rules are satisfied. The advantages 
of this scheme are that it is simple and that the inten- 
sity of the controlling laser pulse is relatively small. The 
main disadvantages are that the efficiency of the scheme 
is sensitive to the pulse area and that the resonance con- 
ditions should be satisfied accurately, especially in the 
case of multiphoton excitation. 

The schemes based on adiabatic rapid passage (ARP) 
rely on adiabatic properties 0, ■ I'^ this case a sin- 
gle chirped laser pulse is used to establish an adiabatic 
regime so that the complete population transfer occurs as 
a result of system evolution along the initially populated 
adiabatic state. The population transfer is not sensitive 
to the pulse area once an adiabatic regime is established, 
so the method is quite robust. However, the pulse area 
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must be much larger than in the case of 7r-pulses. This 
becomes critical if the competing ionization or dissocia- 
tion processes are significant. 

Stimulated Raman adiabatic passage (STIRAP) 
emerged as a very efficient method to provide almost 
complete population transfer in a three-level system |12| | . 
In its simplest form, the scheme involves a two-photon 
Raman process, in which an interaction with a first 
(pump) pulse links the initial state with an intermedi- 
ate state, which in turn interacts via a second (Stokes) 
pulse with a final target state. An important advantage 
is that almost no population is placed into the intermedi- 
ate state, and thus the process is insensitive to any pos- 
sible decay from that state. In more complicated cases 
the method is used to create a maximal coherent super- 
position of several states and complete population 
transfer via several intermediate states The analog 
of the STIRAP which relies on adiabatic properties is 
called chirped adiabatic passage by two-photon absorp- 
tion (CAPTA) 151 

Recently an effective scheme based on the idea of pe- 
riodic chirping was introduced [H [13. Within this 
scheme, a sequence of chirped laser pulses or a laser with 
periodic chirping is used to establish multiple crossings 
between dressed initial and target levels. Effective laser 
control is performed by manipulating the parameters of 
these crossings and adiabatic phase differences between 
the two crossings directly. In its simplest formulation the 
complete population transfer can be achieved by using a 
single quadratically chirped pulse. Because of interfer- 
ence between the crossings, the necessary pulse area of 
each pulse is much smaller compared to the ARP; how- 
ever the method is still quite robust. These advantages 
make the method useful in the field of wave packet con- 
trol 

PI- 

However, the area of the pulse is still larger 
than in the case of control by non-chirped pulses. 

The methods from the second group are based on max- 
imization of the desired output by a step-by-step adjust- 
ment of the controlling field. The main advantage of 
these methods is that controlling field search is performed 
explicitly, so that complicated systems can be controlled. 
These methods can be arranged according to the measure 
of the controlling parameter space. 

In the simplest case, the number of controlling laser 
pulses and their envelopes are fixed. Each pulse is de- 
scribed by a few parameters: peak intensity, central fre- 
quency, chirp rate, central time and full width at half 
maximum (FWHM) . Control is performed by direct opti- 
mization of these parameters. By solving the Schrodinger 
equation, the final population of the desired state can be 
determined as a function of the parameters of the con- 
trolling pulses. Thus, well established numerical proce- 
dures for maximization of this function of many variables 
can be used |2^. It is easy to demonstrate that all the 
schemes discussed above can be reproduced by direct op- 
timization as particular cases. The approach can be con- 
sidered as an optimization with respect to a finite set of 
parameters, i.e., in a multi-dimensional space. 



To design the controlling field without any restrictions 
on its shape, the optimal control theory (OCT) was devel- 
oped milflllllilllll. It is based on the idea that the 
controlling laser pulse should maximize a certain func- 
tional. The basic variational procedure leads to a set of 
equations for the optimal laser field, which include two 
Schrodinger equations to describe the dynamics starting 
from the initial and target state wave functions. The 
optimal laser field is given by the imaginary part of the 
correlation function of these two wave functions. This 
system of equations of optimal control must be solved 
iteratively in general. The reader can find a comprehen- 
sive review of OCT in The main disadvantage of 
this method is that in many cases the generated optimal 
field is hardly possible to realize in an experiment. The 
approach can be considered as optimization with respect 
to a continious set of parameters, i.e., in the fimctional 
space. 

In the case of selective excitation of metastable atomic 
states by laser pulses, a multiphoton interaction is the 
basic mechanism of the process. Such effects as direct 
ionization ^25. .26] , resonant transitions via intermediate 
discrete and continuum stat es [261 l27l 123 . |2^ Isof , above 
threshold ionization |2^ Is^. I33l l34ll , an d electron rescat- 
tering by the atomic core |35l l36l 13^ l3^ are essential. 
An accurate treatment of all these effects requires the 
exact solution of the quantum equations to describe the 
dynamics of the system. Unfortunately, it is impossible 
to introduce some reliable and solvable model equations 
to estimate the parameters of the controlling laser field 
analytically. Thus, the methods of maximization of de- 
sired output should be used. 

In the present work selective excitations of the 
metastable 2s state of hydrogen and the singlet 
metastable ls2s ^S state of helium are studied using both 
the direct optimization of laser parameters and the op- 
timal control theory. To integrate the time-dependent 
Schrodinger equation, the well established close-coupling 
approach based on the pro pert ies of time-dependent in- 
tegral equations |^ |29|, [33, Hj is used. This approach is 
the most suitable for this kind of spectroscopic calcula- 
tion. As electron-atom collisions rapidly destroy atomic 
metastable states, a search of the optimal laser param- 
eters is performed under the requirement that not only 
high final population of desired metastable state should 
be achieved, but also ionization probability should be 
small jL,5]. 

The present paper is organized as follows. In the next 
section the theoretical approaches are summarized, and 
relations between the results of the optimal control cal- 
culations and the recently known schemes are demon- 
strated. In Sec. mil direct optimization of the parame- 
ters of one- and two laser pulses with Gaussian envelope 
is used to introduce a number of simple schemes for se- 
lective excitation of II(2s) and IIe(ls2s ^S). In Sec. IIVI 
the optimal control theory is used to look for other possi- 
ble controlling laser fields. A new mechanism of selective 
excitation by using attosecond pulses is introduced. Sec. 
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contains a summary. 



II. THEORY 

A. Close-coupling approach 

The time-dependent Schrodinger equation for an atom 
in a laser field is written as 



^--Ho + F{t)d{r) 



(1) 



where Hq is the Hamiltonian of the unperturbed atom, 
d(r) is the dipole moment and F{t) is a laser field. 
Atomic units are used throughout the paper, unless oth- 
erwise noted. 

We employ the close-coupling (CC) method on the ba- 
sis of the orthogonal and normalized unperturbed atomic 
wave functions of the discrete and continuum states and 
expand the total wavefunction as 

where v and /i are the indices that represent the inte- 
ger quantum numbers of discrete and continuum states 
respectively. The functions \il}v) and \i}p,E) are the cor- 
responding stationary state wavefunctions, Ejj are the 
energies of the discrete states, E stands for the energy 
of the electron continuum and a^{t) and bf^Eit) are un- 
known coefficients. 

The Hermitian system of coupled equations for the co- 
efficients of the discrete and continuum states follows 
from substituting the expansion Eq. (0) into the time- 
dependent Schrodinger equation Eq. (^) as 

at ^-^ 



V 

-t-F(t) E / U^M' (^^: E') e<''-'''>b^E'{t)dE'. (4) 



The matrix elements U^i,', \J^fj,{E) and U^^'(£^,£") are 
integrals over r-space taken with the atomic dipole mo- 
ment operator. 

The first and the second sums in Eq. Q describe the 
bound-bound and free-bound transitions, respectively. 



These transitions result in a significant redistribution of 
the population of discrete states. This in turn strongly af- 
fects the ionization process 1^ • The first sum 
in Eq. Q describes ionization from all the discrete states 
and the integral term (free-free transitions) describes 
a multiphoton inverse brcmsstrahlung process (within 
quantum mechanical considerations) or rescattering pro- 
cesses (within quasiclassical considerations). These pro- 
cesses play an important role in formation of the pho - 
toelectronic spectra at high energies HE IHE l38| . 
However, since the free-free matrix elements practically 
do not affect the discrete state amplitudes a^, (t) , it is 
possible to neglect them for the evaluation of {t). 
The transitions neglected in this approximation are the 
third-order (bound- free- free-bound) ones. The role of the 
free-free transitions was carefully investigated in Refs. 

and this assumption was confirmed 

to work well. 

The close-coupling equations Eqs. (|3I4|I with the free- 
free transitions neglected can be solved by discretizing 
the continuum or by employing the recently developed 
approach based on the properties of the time-dependent 
integral equations without any discretization of the con- 
tinuum [23, 12E HE 1^ • As the nonresonant interaction 
between the lower and highly excited discrete states is 
small, it is possible to adjust the number of discrete 
states involved in expansion Eq. (01 so that the results of 
the calculations do not change. Although, strictly speak- 
ing, this procedure does not prove the convergence of the 
close-coupling approach, it is still widely used in coUi- 
sional physics [y, |33 • 

In the present study we consider excitation and ion- 
ization of a single atomic electron in the field of linearly 
polarized laser pulses. 



B. Methods of direct optimization of laser 
parameters 

In the simplest case the controlling field is assumed to 
be a sequence of laser pulses and is written as 



N 



i=i 



X sm 



(5) 



where / (t, At) is a fixed envelope function with FWHM 
equal to r and centered to achieve maximum at time 
difference /S.t — Q. In the present study we focus on the 
laser pulses Eq. Q with a Gaussian shape: 



/ (r, At) = exp 



4 In 2 



(6) 



The system is controlled by changing peak amplitudes 
Fj, frequencies Wj, chirp rates a^, central times tj, and 
FWHM's Tj of the component pulses. By solving the 
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Schrodinger equation one can determine the final popula- 
tion of a selected state as a function of these parameters. 
To find the maximum of this function, the conjugate gra- 
dient search method is used here. The gradient of 
the final population of the target state with respect to 
the laser parameters is calculated numerically using the 
finite-difference approximation. The optimization starts 
from some initial guess parameters. 

As different laser pulse parameters have different di- 
mensions, the parameters should be put on a common 
ground by using dimensionless units, or the optimizations 
with respect to different types of parameters should be 
performed separately. In the first case the convergence of 
the optimization procedure could depend on the dimen- 
sionless units used. Indeed, the efficiency of the selec- 
tive excitation by non-chirped pulse strongly depends on 
whether or not the resonant conditions are achieved, so 
that by using common atomic units one should find the 
maximum of the function that strongly depends on one 
group of arguments (frequencies) and only weakly on the 
other arguments (intensities, durations etc.). This exam- 
ple demonstrates that proper dimensionless units should 
be used to ensure reasonable convergence. Moreover, the 
derived optimal parameters could depend on how that 
units are defined. 

In the present study the optimizations with respect to 
different types of parameters are performed separately so 
that each step of the optimization procedure should con- 
tain a set of optimizations with respect to the parameters 
for all pulses that have the same dimension. Generally, 
some laser parameters should be linked with each other to 
achieve a maximum of the desired output. If the chirping 
rate is zero, the optimal intensity and duration of each 
pulse are related to each other by some pulse area con- 
servation rule. To avoid this uncertainty, we assume the 
durations of the pulses to be equal. After the optimal 
parameters of the controlling sequence of laser pulses are 
found, we can adjust intensities and durations to fit the 
specifics of the experimental technique. 

The order of optimizations performed at each step of 
the direct search affects the efficiency, convergence and 
resulting controlling scheme. However, since the laser 
field, Eq. 0, has only five different types of parame- 
ters, the number of possible nonequivalent orders of op- 
timization is limited. Table presents the correspon- 
dence between the controlling schemes generated by the 
direct optimization method with different orders of opti- 
mizations performed at each step and previously known 
schemes. Since the last schemes have been established 
without bound-free transitions taken into account the 
continuum states are not included in the test calculations 
of Tabled The continuum states are, however, taken into 
account at all other calculations, reported in the present 
paper. 

The durations of all the pulses are assumed to be equal 
and do not change during optimization. Initial param- 
eters used in all the calculations are 1 GW/cm^ for in- 
tensities of all pulses, frequencies are estimated from the 



TABLE I: Controlling schemes generated by the direct op- 
timization method with different orders of optimizations per- 
formed at each step. 
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Case 


— 

Number 




Order of 


Known 


No. 


of pulses 


optimizations 


particular cases 


1 


1 


1. 
2. 


Intensities 
Frequencies 


TT- pulses 


2 


1 


1. 
2. 


Intensities 
Chirping rates 


ARP 


3 


2 


1. 
2. 


Central times 
Intensities 


Direct excitation 
or STIRAP 






3. 


Frequencies 


4 


2 


1. 
2. 


Central times 
Intensities 


Direct excitation 
or CARTA 






3. 


Chirping rates 


5 


Several 


1. 
2. 
3. 


Frequencies 
Central times 
Intensities 


Periodic sweeping 
of laser intensity 


6 


Several 


1. 
2. 
3. 


Chirping rates 
Central times 
Intensities 


Periodic sweeping 
of laser frequency 



data for the unperturbed atom, chirping rates are zero 
for all the pulses, centers of the pulses coincide in the 
cases 1-4 and are equally distant with the 2-FWHM shift 
in cases 5 and 6. Our calculations show that number of 
steps required to achieve convergence is about one and a 
half the number of laser pulses used in the scheme. 

The procedure of direct optimization of laser parame- 
ters can also be used to calculate intrinsic parameters of 
the atomic system. For example, by fixing the intensity 
of the laser and calculating the optimal frequency of se- 
lective excitation of different target states, the dynamic 
Stark shifts of these states can be easily estimated. 



C. Methods of the optimal control theory 

The dynamics of the atom in the laser field is described 
by the time-dependent Schrodinger equation The 
initial state wave function \^i{t)) is specified at time 
t = 0. The goal of control is to design such an exter- 
nal field F(t) that the wave packet calculated with Eq. 

up to time t = T is close enough to the desired target 
state wave packet |$f(t)). One of the most natural and 
fiexible approaches to design such a field is the optimal 
control theory . It is based on the idea that the con- 
trolling laser pulse should maximize a certain functional. 
The procedure leads to a set of equations for the optimal 
laser field, which include two Schrodinger equations to 
describe the dynamics starting from the initial and tar- 
get state wave packets. The optimal laser field is given 
by the imaginary part of the correlation function of these 
two wave packets. This system of equations of optimal 
control must be solved iteratively in general starting from 
some initial guess field. 

A number of algorithms to realize this idea have been 
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developed [Tl[^El!25,'2?|. One of the most effective is 
the algorithm by Zhu, Botina and Rabits |23i| (ZBR algo- 
rithm), which is developed to solve the following system 
of the optimal control equations: 



i^-i7o + r(i)d(r) 



|0(t))=O, 10(0)) = !$,), (7) 



\x{t))=0, |x(T)) = |$*), (8) 



pulse for two-photon excitation of hydrogen (wavelength 
is around 240 nm) and helium (wavelength is around 120 
nm). The frequency of the laser was optimized for each 
pair of the arguments to compensate level shifts due to 
the dynamic Stark effect and to maximize the output. 

For the case of selective excitation of H(2s) the max- 
imal values of the target state population achieved are 
18 % to 20 %. Unfortunately, the corresponding ioniza- 
tion probabilities are large, 27 % to 30 %. The maxi- 
mal differences between the target state population and 
ionization probability are found for the following set of 
parameters: 



F{t) ^ — Im mt) \x{t)) (x(i)l d(r) |0(t))] . (9) 
7 

Here 7 is a positive parameter chosen to weight the sig- 
nificance of the laser energy (2^]. To integrate Eq. (0 
forward in time and Eq. (jSJ backward in time, we employ 
the close-coupling method discussed in subsection III Al 

To launch the ZBR algorithm, an initial field should 
be specified. In molecular dynamics control, which is 
the main area of the application for the ZBR algorithm, 
a zero initial field is sufhcient for most cases. However, 
since in the present study initial and target states are sta- 
tionary, a zero initial field cannot be used. Instead we use 
a single laser pulse of Gaussian shape as an initial field 
and generate different optimal control fields by changing 
its parameters. Our calculations show that three to five 
iterations of the ZBR algorithm are sufficient to achieve 
the convergence. 



III. SELECTIVE EXCITATION OF THE 
METASTABLE ATOMIC STATES BY 
FEMTOSECOND PULSES 

In the present section selective excitations of H(2s) and 
He(ls2s ^S) by one and two femtosecond laser pulses of 
Gaussian shape are studied. The case of chirp pulses 
miSini or periodic chirping [H [13 is not considered 
as the pulse area required is larger than in the case of 
non-chirped pulses. This in turn significantly increases 
the ionization. However, in some cases when the ioniza- 
tion process is suppressed, periodic chirping can be used 
to improve the robustness of the scheme. These will be 
discussed in a future publication. 



A. Selective excitation of H(2s) and He(ls2s ^S) by 
a single laser pulse 

In the simplest case a two-photon excitation of H(2s) 
and He(ls2s ^S) can be achieved by using a single laser 
pulse. The process can be affected by changing only three 
parameters of the laser: frequency, duration, and inten- 
sity. To study this case we calculate the final populations 
of the target metastable states and ionization probabil- 
ities as functions of intensity and duration of the laser 



FWHM [fs] = 2.25 x 10^'^/ 1 [W/cm^], (10a) 

LU [eV] = 5.1021 eV 

+ 2.34 X 10"^^ X / [W/cin] (10b) 

where / is the pulse intensity. The target state popula- 
tion and ionization probability are found to be 8.4 % and 
4.2 %, respectively, for this set of FWHM and ut. Equa- 
tions (|10|l are obtained by a fit performed for the range 
of intensities from 1 TW/cm^ to 20 TW/cm^. 

The efficiency of the process is low because of sig- 
nificance of the one-photon ionization from the target 
metastable state. As the photon energy is about 5.10 
eV to 5.14 eV (depending on pulse intensity), this pro- 
cess populates the continuum states with low energies 
of about 1.7 eV, so that the corresponding bound-free 
transition matrix elements are large. 

An opposing situation is observed in the case of selec- 
tive excitation of He(ls2s ^S). The best values of 65 % 
for the target state population and 26 % for ionization 
probability are found for the following parameters: 

FWHM [fs] = 7.12 X 10^^/ [W/cm^], (11a) 

uj [eV] = 10.3075 eV 

-I- 1.54 X 10"^^ X / [W/cin]. (lib) 

Eqations Hll|l are obtained by a fit performed for the 
range of intensities from 50 TW/cm^ to 400 TW/cm^. 

The photon energy is about two times higher then in 
the case of hydrogen (10.3 eV to 10.4 eV, depending on 
pulse intensity). The process is effective. One-photon 
ionization from the target metastable state populates 
the continuum states with energies of about 6.3 eV and 
corresponding bound-free transition matrix elements are 
small, so that ionization does not undermine the process. 



B. Selective excitation of H(2S) by two laser pulses 

To improve the efficiency of selective excitation of the 
metastable state of hydrogen, we increase the number 
of controlling parameters by introducing two controlling 
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-■— 40fs a 40fs 
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-•— 80fs o aofs 

lOOfs O lOOfs 




Frequency of the "blue" pulse [eV] 

FIG. 1: Selective excitation of H(2S) by two laser pulses; 
(top) final populations of 2s metastable state and ionization 
probability, (middle) time shift between the "blue" and "red" 
pulses and (bottom) intensities of the pulses given as functions 
of frequency of the "blue" pulse. 



laser pulses with different frequencies instead of a sin- 
gle pulse. For simplicity, we call the pulse with lower 
frequency as "red" pulse and the pulse with higher fre- 
quency as "blue" pulse. Frequencies of the blue and the 
red pulses should be linked to fit the energy difference be- 
tween the ground and the metastable states. As optimal 
intensity and duration of each pulse are approximately 
related to each other by a pulse envelope area conser- 
vation rule, in the present study we fix durations of the 
pulses to be equal. 

Figurenpresents the final population of the metastable 



state and ionization probability of hydrogen as functions 
of frequency of the blue pulse for different durations of 
the pulses with the frequency of the red pulse, intensities 
and central times of both pulses adjusted to maximize 
the population of the metastable state. 

If the frequency of the red pulse becomes lower than 
the threshold of the one-photon ionization of the 2s 
metastable state, the efficiency of the selective excitation 
grows dramatically. This is caused by an abrupt decrease 
of ionization rate, which makes it possible to apply more 
intense laser pulses and achieve higher selective excita- 
tion without increasing the ionization. Unfortunately, 
for hydrogen the one-photon ionization threshold for the 
2s state is equal to one third of the energy difference be- 
tween the Is and 2s states, and the simple case, when 
the blue pulse is produced from the single red pulse by 
the second harmonic generation, is still not good enough. 
Acceptable results, however, can be obtained for the fre- 
quencies of the blue pulse slightly higher and frequencies 
of the red pulse slightly lower than the threshold val- 
ues. An interesting observation comes from the analysis 
of the time shift between the two optimized controlling 
pulses. If the frequencies of the pulses are far from being 
in resonance with some bound state, their centers coin- 
cide. However, if the frequency of the blue pulse is close 
to the frequency of 3p — > Is transition, the centers of the 
blue and the red pulses separate. In this case the direct 
two-photon excitation process transforms into STIRAP 
[3 (Table n case 3). 

Consider three examples in more detail: 

Example 1. Final populations: 2s - 59 %, continuum 
- 21 %. Parameters of the pulses. Freqencies: blue - 
7.12786 eV, red - 3.13629 eV; Intensities: blue - 5.8 x 
10^2 iy/cm2, red - 1.2 x lO^^ W/cm^; FWHMs are 40 fs 
for both pulses; Pulse centers coincide. The frequency of 
the red pulse is below but close to the threshold of the 
one-photon ionization of the H(2s) state. Although the 
frequencies are adjusted to avoid resonances with highly 
excited discrete states, 7p and 8p states are populated at 
the level of 9 % and 3 % respectively. 

Example 2. Final populations: 2s - 92 %, con- 
tinuum - 7 %. Parameters of the pulses. Freqencies: 
blue - 9.09943 eV, red - 1.09933 eV; Intensities: blue - 
4.8 X 10^2 W/cm'^, red - 1.6 x lO^^ W/cm"^; FWHMs: 
blue - 40 fs, red - 100 fs; Pulse centers coincide. In this 
case durations of the pulses are different and efficiency 
of the process is less sensitive to the shifts between the 
pulse centers. 

Example 3. Final populations: 2s - 91 %, con- 
tinuum - 8 %. Parameters of the pulses. Freqencies: 
blue - 12.0878 eV, red - 1.88656 eV; Intensities: blue - 
3.7 X 10^2 w/cm'^, red - 4.0 x 10" W/cm'^; FWHMs are 
100 fs for both pulses; blue pulse is delayed for 58 fs. The 
frequencies of the pulses are close to be in resonance with 
3p state. In this case the red pulse comes before the blue 
pulse. This represents the STIRAP process [T^ . 
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FIG. 2: The optimal control field generated for selective exci- 
tation of the ls2s state of helium (top), its Fourier trans- 
form (middle) and time variation of the metastable state pop- 
ulation and ionization probability of the controlled atom (bot- 
tom). Top and middle: dashed line is the OCT result, solid 
line is the fitted sequence of attosecond pulses. Insert: initial 
guess field. Bottom: solid line is the population of the ls2s 
state, dashed line is the ionization probability. 



IV. SELECTIVE EXCITATION OF THE 
METASTABLE ATOMIC STATES BY 
ATTOSECOND PULSES 



In the previous sections the controlling laser field was 
assumed to be composed of a few femtosecond laser pulses 
with Gaussian envelopes. To study other possibilities 
to achieve selective excitation we use the optimal con- 
trol theory, as described in Sec. Ill CI As the controlling 
schemes introduced in the previous section for long du- 
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FIG. 3: Final populations of the metastable states and ion- 
ization probabilities obtained with the optimal field generated 
from a single Gaussian laser pulse given as a functions of the 
parameters of initial guess pulse. Top line: Results obtained 
for selective excitation of H(2s) with the frequency of initial 
guess pulse is 10.2 eV. Bottom line: Results obtained for se- 
lective excitation of H(ls2s ^S) with the frequency of initial 
guess pulse is 10.3 eV. 



rations can hardly be improved, we focus on the case of 
short excitation times. 

A typical example is shown in Fig. |21 It presents the 
optimal control field generated by the ZBR algorithm for 
selective excitation of ls2s state of helium (top) and 
the time variation of population of the target metastable 
state and ionization probability during the interaction 
(bottom). The calculations are performed starting from 
the Gaussian pulse of intensity 80 TW/cm^, frequency 
10.3 eV and FWHM 2.5 fs (full duration is 5 fs), which 
is shown in the insert of Fig. [3 The weight parameter 
is taken to be 7 = 0.1. The final population of the tar- 
get state is 77 % with ionization probability 7 %. The 
convergence is achieved in four iterations. 

During optimization the laser field completely trans- 
forms into a sequence of attosecond pulses. The main 
perturbations finish after 2 fs, so the duration of active 
control is also changed. Contrary to the case of con- 
trolling femtosecond pulses, the state populations now 
change abruptly. Physically, the sequence of attosecond 
laser pulses affects the atomic electrons as a series of 
pushes. Each push displaces the electrons for a small 
distance so they can not move far from the atomic core, 
and this suppresses ionization. 

To simplify the optimal control field obtained with the 
ZBR algorithm we approximate it by a sequence of Gaus- 
sians: 



N 



A, 



F(t)=^^expM(t-i,) 



(12) 



If the duration of each pulse is short enough, the momen- 
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turn supplied to the electrons by each attosecond pulse 
is given by the pulse area Aj and does not depend on 
(3. To calculate the parameters Aj and tj we roughly ap- 
proximate the controlling field by sequence H12|l and use 
the method described in Sec. Ill BI to optimize Aj and tj. 
The resulting sequence of attosecond pulses is presented 
in Fig. 121 (top) as a bold line. The calculation is per- 
formed with the width parameter /3 = 0.04 (FWHM = 
20 as). The final population of the target state and the 
ionization probability for that field are 81 % and 14 %, 
respectively. 

Figure |3| presents final populations of the metastable 
states and ionization probabilities obtained with the op- 
timal fields generated from femtosecond pulses of differ- 
ent intensities and durations. Each point of the plot 
represents an effective controlling sequence of attosec- 
ond pulses. One can see that the results depend on the 
parameters of the initial guess pulse in a very compli- 
cated way. The parameters for some effective controlling 
sequences of attosecond pulses Eq. (IT^ calculated by 
optimal fitting of the parameters Aj and tj are presented 
in the Table |ni 

Recently, q number of schemes to generate trains of 
attosecond laser pulses have been proposed and realized 
experimentally (see review 113). The basic idea is the 
production of a comb of equidistant frequencies in the 
spectral domain with controlled relative phases. As a re- 
sult trains of x-ray laser pulses of attosecond duration are 
obtained. The shapes and carrier frequency can be mod- 
ified by eliminating certain harmonics. Also, waveforms 
containing optical bursts approaching one cycle have 
been successfully synthesized by superposing five phase- 
controlled sidebands ^Tj- Trains of non- modulated at- 
tosecond pulses could be also generated b y la ser-plasma 
interaction in the relativistic regime (see i42i| and refer- 
ences therein.) 

To generate sequence of pulses close to ones shown in 
Fig. 2, the spectral components with the frequencies 7.33 
eV, 22.0 w 3 X 7.3 eV, 29.3 w 4 x 7.3 eV and 44.0 w 6 x 7.3 
eV should be syncronized. Four main harmonics taken 
from equidistant spectral modes should be superposed. 
Thus, for principal frequency of 7.33 eV one should take 
1st, 3rd, 4th and 6th harmonics; for principal frequency 
of 3.66 eV 2nd, 6th, 8th and 12th harmonics, respectively, 
and so on. 



V. CONCLUSIONS 

In the present study the controlling field search is per- 
formed using a well established close-coupling approach 
for solving the time-dependent Schrodinger equation, so 
that all essential multiphoton effects are treated accu- 
rately. It is demonstrated that optimization procedures 
provide an effective technique to design laser fields for se- 
lective excitation of metastable atomic states. Not only 
do they reproduce the recently developed schemes of laser 
control as particular cases but also they introduce new 
ones. 

An efficient selective excitation of II(2s) and He(ls2s 
^S) can be achieved by using one and two femtosecond 
laser pulses. Frequencies of the pulses should be ad- 
justed properly to suppress single-photon ionization from 
the target metastable state. Optimal populations of the 
atomic states and the corresponding parameters of the 
laser pulses are calculated as functions of preferable fre- 
quencies and durations of the pulses. These make it pos- 
sible to choose the parameters of the controlling laser 
pulses that fit the capabilities of the present experimen- 
tal techniques. 

A new way to achieve selective excitation of metastable 
atomic states by using sequences of attosecond pulses is 
introduced. This is important because of recent progress 
in production of attosecond pulses. While in the present 
study the parameters of the pulses are calculated us- 
ing the optimal control theory, the direct optimization 
of high harmonics can be used to control the process in 
the future. 
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